a = 0.1;
b = 0.3;
r = 1.0;


N = 1000;

theta = linspace(0,2*pi,N);
x = a + r*cos(theta);
y = b + r*sin(theta);
z = x.*y;

l1norm = abs(x) + abs(y) + abs(z);

atan2(1,1);

solution_theta = atan2(u(2)-b,u(1)-a);
while solution_theta < 0
    solution_theta = solution_theta + 2*pi;
end

plot(theta,l1norm)
hold on
scatter(solution_theta, results(k).u_l1norm, 20, 'r', 'filled');